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This  paper  deals  with  Fermi  systems  at  zero  temperature  whose 

unperturbed  density,  n  (r),  varies  slowly  on  a  scale  Z  »  p_^,  where 

O'*  r 

pp  is  a  typical  local  Fermi  momentum.  We  wish  to  calculate  the 
linear  density  response  of  such  a  system  to  an  external  scalar 
potential,  v®(r,t),  which  we  write  as 

Vje(r,t)  -  — — ,  /  v*(q,u)ei^S’£  "  Wt^dq  du>  (1) 

(2TT)4 

For  the  time  being  we  consider  only  neutral  systems  with  short  range 
forces;  extension  to  charged  systems  poses  no  significant  problems. 

For  several  (overlapping)  regimes  of  q  and  w  the  necessary  theory 
is  already  available.  In  view  of  the  slow  spatial  variation  of  nQ(r), 
a  plausible  ansatz  for  the  density  response  is  the  local  density 
approximation 

ni  (£,<*>)  *  /  X  (£-£*,  w;  nQ(r))  v,e(r',  w)  dr'  ,  (2) 

where  x(r"r,«  w»  nQ)  is  the  dynamic  response  function  of  a  system 

of  uniform  density,  n  .  This  ansatz  is  in  fact  valid  under  the 

o 

following  circumstances: 

(I)  w  «  Vpq  (Vp  •  Fermi  velocity).  In  this  case  X  is  to  a 
good  approximation  the  static  x>  which  is  known  to  have  a  short  range 

of  the  order  of  p”1  or  of  the  range  of  the  forces  (region  (I)  in  Fig.  1) 

(II)  w  »  Vp/i.  Again  x  bas  a  short  range  of  the  order  of  Vp/u>  or 

Pp1  (region  (II)  in  Fig.  1). 

-1  e 

(III)  q  »  Z  .  Here  q  denotes  a  characteristic  wave-number  of  Vj 

The  effective  range  of  x  is  no  larger  than  q  1  (region  (III)  in  Fig.  1 jV 


2 


It  can  be  seen  from  Fig.  1,  that  there  is  a  region,  R,  shown 
enclosed  by  heavy  lines,  which  is  not  covered  by  the  above  three 
regions.  Here  we  develop  a  theory  which  is  valid  in  the  following  regime, 

(IV)  q  «  pF,  u)  «  pF  (region  (IV)  in  Fig.  1,  shown  enclosed  by 

dashed  lines).  This  is  the  region  of  validity  of  the  classical 

(2  ) 

Landau  theory  of  Fermi  liquids.  If  £  is  large  enough,  region  (IV) 
includes  R,  and  therefore  our  theory  compliments  the  local  density 
approximation  (2). 

Our  theory  is  the  generalization  of  Landau  theory  for  systems 
of  slowly  varying  density.  Since  the  Landau  theory  for  homogeneous 
systems  is  valid  for  long  wavelength  perturbation,  it  is  plausible 
that  it  can  be  extended  to  inhomogeneous  systems  whose  unperturbed 
density  varies  on  a  scale  £  large  compared  to  p~^. 

For  a  homogeneous  system,  subject  to  an  external  potential 
Vje(r,t),  the  Landau  transport  equation  for  the  quasiparticle  dis¬ 
tribution  function,  n.(p,r,t),  is  obtained  from  the  classical  Hamilton's 
equations  using  the  effective  Hamiltonian 

H(p,r,t)  -  e(p)  +  v,*(r,t)  +  f  ,  n^p'.r.t)  .  (3) 

-  -  -  -  _i  PP  1  -  - 

P 


Spin  indices  have  been  suppressed.  Here 


2m*  2m* 


(A) 


is  the  one  particle  energy  and  the  last  term  arises  from  short  range 
quasi-particle  Interactions.  Our  generalization  for  an  inhomogeneous 
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system  is 

H(p,r,t)  -  e(p,r)  +  Vj®(r,t)  +  X  f  ,  (r)  n  (p*,r,t)  (5) 

~  ~  pi  PP  ^  1  ^  ~ 

«%* 

where  , 

2  Pp  <E> 

e(p.r)  — 2 - - -  ;  (6) 

2m* (r)  2m* (r) 

f  ,  (r),  m*(r)  and  p_,(r)  are  equal  to  the  homogeneous  Fermi  liquid 
PP  ~  *  ~ 

parameters  f  , ,  m*  and  p_  evaluated  for  the  local  density,  n  (r). 

pp  r  o  «. 

*>**»  /  » 

This  ansatz  is  in  the  spirit  of  local  density  functional  theory  ' 

where  all  local  physical  quantities  are  determined  by  the  local  density. 

The  single  particle  energy  e(p,r)  can  be  derived  from  the  theory 

(||) 

of  single  particle  properties  of  the  inhomogeneous  electron  gas  . 

It  can  also  be  understood  as  follows.  The  local  wave  number  of  a 
quasiparticle  of  vanishing  excitation  energy  is  equal  to 


(This  relationship  is  unaffected  by  interactions).  In  a  system  in 
equilibrium,  quasiparticles  of  vanishing  excitation  energy  must  have 
the  same  single  particle  energy  at  any  two  points  r}  and  r2.  In  view  of 
(7)  this  is  correctly  described  by  the  expression  (6).  Further,  at  a 
given  r,  the  momentum  dependence  of  c  for  small  p  -  p_(r)  must  be 
given  by  the  local  value  of  m*(r)  as  is  the  case  for  our  expression 
(6). 

The  external  perturbation,  v}e  and  the  interaction  term  of  H 
are  the  same  as  in  the  homogeneous  Landau  theory,  except  that. 
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because  of  the  short  range  of  the  interaction,  f  , (r),  is  not  a 

PP  - 

constant  but  corresponds  to  the  local  density. 

The  collisionless  transport  equation  is 

+  p  •  Vp  n  +  r  •  7r  n  -  °  (8) 

where  p  and  r  are  determined  by  Hamilton's  equations, 

p  -  -  V  H  ;  r  -  V  H  (9) 

~  r  p 

We  write 

n(p,r,t)  -  n  (p,r)  +  n  (p,r,t)  ,  (10) 

O  %  %  * 

substitute  (9)  into  (8)  and  linearize:  . 


This  integro-dif ferential  equation  can  be  solved  by  a  variety  of 
techniques.  Here  we  regard  it  as  describing  the  flow  of  classical 
particles  under  the  action  of  the  Hamiltonian  e(p,r)  and  with  a  (self- 
consistent)  source  S, (p,r,t)  equal  to  the  right  hand  side.  Consequently 
we  are  led  to  define  the  following  Green's  function 

G(p,r  ;  p'r'  ;  t-t')  -  5(p'  -  P  (p,r,t')) 

6(r'  -  Q  (p,r,t'))  0(t-t') 

%  %  %  % 


(12) 
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where  Q,P  describe  the  trajectory  of  a  particle  moving  under  the 
action  of  e(P,Q),  which  was  created  at  time  t'  at  r'  with  momentum 
p'  and  arrives  at  r  with  momentum  p  at  time  t.  The  trajectories  are 
given  by 

P  -  -  VQ  e(P,Q),  Q  -  Vp  e(P,Q)  (13a) 

and 

P(t)  “  P  ,  Q(t)  -  r  (13b) 

The  formal  solution  of  Eq.  (11)  is  then  given  by 


n^p.r.t)  -  nj(p,r,t)  +  /G(p,r;p' ,r' ;t-t ' )SJ(p' ,r' ,t ' )dp'dr*  dt’ 

(14) 

where  n}  is  a  solution  of  the  homogeneous  equation.  If  we  now  consider 
a  single  time  Fourier  component  of  v®,  and  the  corresponding  density 


change 


n  -  o  . 


then  the  equation  (14)  can  be  reduced  to 


n1(P*£»w)  *  (p,r,«)  -iw/Sj  (p1  ,r’  ,u>)  G(p,r;p,r’  ,co)  dp'dr'J 


where 


Sl(p,rfaj)  -  v*(r,u)  +S  f  , (r)n.(p',r  ,u>) 
~  ~  1  -  _i  PP  -  ‘  ~  ~ 


and 
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(22) 


7 


*V 


|| 

here  X  is  the  response  function  of  a  homogeneous  system  of  density 

LH  LI 

oq;  x  is  the  same  quantity  calculated  from  Landau  theory;  and  \ 

is  the  response  function  calculated  in  the  present  extension  of  Landau 

theory  to  inhomogeneous  systems.  Provided  £  is  large  enough,  this 

expression  becomes  arbitrarily  accurate  for  all  external  perturbations 

Vj  (r,w)  regardless  of  the  value  of  oj  and  the  characteristic  wave 

numbers  q.  For,  in  regions  (II)  and  (III),  the  last  two  terms  cancel 
|| 

and  X  becomes  an  accurate  representation;  and  in  the  remaining  region 
(q  <  £  *  and  oj  <  Vp£  ^),  the  first  and  last  term  cancel  and  x^*  becomes 
an  accurate  representation. 

The  expression  (22)  defines  the  response  function  as  a  functional 
H 

of  the  density:  X  is  determined  by  nQ  and  theoretical  approximations 
to  it  are  available^;  XLH  is  well  known  in  terms  of  the  density- 
dependent  Landau  parameters^;  and  x^*  can  be  calculated  with  the 
same  Landau  parameters  by  our  theory.  Thus  the  expression  (22)  can 
be  regarded  as  a  dynamic  response  function  in  density  functional  theory. 

For  a  non-interacting  one-dimensional  system,  Eq.  (16)  can  be 
solved  and  leads  to  the  following  result  for  x: 


X(x,x',u) 


1  5(x-x’)  ioi _  iO)t(x,x’) 

PF(x)  PF(x)pp(x')  6 


(23) 


where  t(x,x’)  is  the  time  of  transit  between  positions  x  and  x'.  This 
formula  has  been  checked  by  an  Independent  calculation  using  the  BWK 
approximation.  It  is  interesting  to  note  that  for  u>  *  0,  this  X  is 
strictly  local;  for  w  «  vFq  the  second  term  is  effectively  small 
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compared  to  the  first  but  of  long  range;  and  for  larger  u  the  second 
term  increases  and  becomes  effectively  of  shorter  range. 

Extensions  of  this  theory  to  other  forms  of  external  perturba¬ 
tion  are  under  study. 

Although  strictly  speaking  our  theory  applies  only  to  systems 
of  very  slow  spatial  variation,  it  should  be  interesting  to  apply 
it  to  real  systems  whose  density  varies  on  an  atomic  scale. 

This  work  has  been  supported  in  part  by  the  National  Science 
Foundation  and  the  Office  of  Naval  Research. 
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FIGURE  CAPTION 

Fig.  1.  Region  R  is  enclosed  by  heavy  lines  and  region  (IV)  by  the  dashed 

lines.  Local  density  theory  is  valid  in  regions  (I),  (II)  and  (III) 
but  fails  in  region  R;  present  theory  is  valid  in  region  (IV).  a  is 
a  numerical  constant,  large  compared  to  unity;  and  8  is  a  constant  small 


compared  to  unity. 


